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Abstract 

Helium and hydrogen recombination lines observed in low-metallicity, extragalactic, 
Hll regions provide the data used to infer the primordial helium mass fraction, Yp. In deriving 
abundances from observations, the correction for unseen neutral helium or hydrogen is usually 
assumed to be absent; i.e., the ionization correction factor is taken to be unity (icf = 1). In a 
previous paper (VGS), we revisited the question of the icf for Hn regions ionized by clusters 
of young, hot, metal-poor stars, confirming earlier work which had demonstrated a "reverse" 
ionization correction: icf < 1. In VGS the icf was calculated using more nearly realistic models of 
inhomogeneous Hn regions, revealing that for those Hn regions ionized by young stars with "hard" 
radiation spectra the icf is reduced even further below unity compared to homogeneous models. 
Based on these results, VGS suggested that the published values of Yp needed to be reduced by 
an amount of order 0.003. As star clusters age, their stellar spectra evolve and so, too, will their 
icfs. Here the evolution of the icf is studied, along with that of two, alternate, measures of the 
"hardness" of the radiation spectrum. The differences between the icf for radiation-bounded 
and matter-bounded models are also explored, along with the effect on the icf of the He/H ratio 
(since He and H compete for some of the same ionizing photons). Particular attention is paid 
to the amount of doubly-ionized helium predicted, leading us to suggest that observations of, or 
bounds to, He ++ may help to discriminate among models of H n regions ionized by starbursts of 
different ages and spectra. We apply our analysis to the Izotov h Thuan (IT) data set utilizing 
the radiation softness parameter, the [OIII]/[OI] ratio, and the presence or absence of He ++ to 
find 0.95 < icf < 0.99. This suggests that the IT estimate of the primordial helium abundance 
should be reduced by AY« 0.006 ± 0.002, from 0.244 ± 0.002 to 0.238 ± 0.003. 
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1. Introduction 

As the second most abundant element in the universe, whose abundance even today is 
dominated by its early universe production during Big Bang Nucleosynthesis (BBN), 4 He plays 
a key role in testing the consistency of the standard, hot, big bang model of cosmology and 
in using primordial nucleosynthesis as a probe of cosmology and particle physics (for a recent 
review and further references, see Olive, Steigman & Walker 2000). To minimize the unavoidable 
contamination from post-BBN production of 4 He, low-metallicity, extragalactic H n regions are 
the targets of choice for deriving the primordial abundance. Observations of large numbers of such 
Hn regions have led to estimates of the primordial helium mass fraction, Yp, whose statistical 
uncertainties are very small, ~ 1% (e.g., see, Olive &; Steigman 1995 (OS), Olive, Skillman, & 
Steigman 1997 (OSS), Izotov, Thuan & Lipovetsky 1994, 1997 (ITL), Izotov & Thuan 1998 (IT), 
Peimbert, Peimbert, & Ruiz 2000 (PPR)). Now, more than ever before, it is crucial to attack 
potential sources of systematic error. In a previous paper (Viegas, Gruenwald, & Steigman 2000 
(VGS)) we explored the size of the ionization correction for unseen neutral hydrogen and/or 
helium in Hn regions ionized by young, hot, metal-poor stars. In VGS models of more nearly 
realistic, inhomogeneous Hn region models were compared with the data set assembled by IT, 
leading to the conclusion that, on average, icf < 1 for the IT Hn regions. Since IT (as well as 
almost all other analyses of Hn region observations) assumed icf = 1, VGS suggested that their 
derived value of Yp needed to be reduced by « 0.003, from Yp « 0.244 to Yp « 0.241. Some 
support for this lower value comes from the recent observations by PPR of a relatively metal-rich 
SMC Hii region for which they derive Y(SMC)= 0.2405 ± 0.0017, as well as from an independent 
theoretical study by Sauer & Jedamzik (2001; hereafter SJ). 

In VGS models were calculated of Hn regions ionized by the radiation from starbursts of 
two different ages (Cid-Fernandes et al. 1992): t = and t = 2.5 Myr. The predictions for these 
two cases are in stark contrast, with icf < 1 for the young star cluster and icf > 1 for the older 
one. The difference is easy to understand as a consequence of the evolving radiation spectrum. 
The young star cluster has many massive, hot stars which produce a "hard" spectrum with many 
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helium-ionizing photons. As the star cluster ages, the more massive stars die, and the spectrum 
softens, reducing the relative number of helium to hydrogen ionizing photons. However, the 
evolution of the radiation spectrum from a cluster of stars is not monotonic. As time goes by, 
stellar winds strip away the atmospheres of lower mass stars exposing their hotter interiors and 
the resulting stellar spectra harden, producing relatively more He and He + ionizing photons. 
Eventually, even these stars end their lives and the spectrum softens, the overall flux of ionizing 
radiation decreases, and the H n region fades away. 

When comparing with the IT data, in order to distinguish between Hn regions ionized by 
"hard" spectra for which icf < 1 and those ionized by "soft" spectra leading to icf > 1, VGS 
utilized the Vilchez-Pagel "radiation softness parameter" (Vilchez & Pagel 1998) 

v = (^l)(^l) (1) 

finding that the IT data prefer a hard spectrum (log i] < 0.4), favoring a "reverse" ionization 
correction (icf < 1). 

In this paper the evolution of model Hn regions is studied tracking the changes in the icf, 
the radiation softness parameter rj, the abundance of doubly ionized helium (He ++ /H + ), and an 
alternate radiation softness parameter proposed by Ballantyne, Ferland, & Martin (2000; hereafter 
BFM), the ratio [OIII]A5007/[OI]A6300. In §2 we follow the evolution of the radiation spectrum 
using the starburst models of Cid-Fernandes et al. (1992), concentrating on the relative changes 
in the hydrogen and helium ionizing fluxes and their implications for the evolution of the icf, 
77, [OIII]/[OI], and the HeII(4686)/H/3 ratio. Since hydrogen and helium compete for the same 
ionizing photons, the sizes of the H° and He zones may depend on the "true" He/H ratio. This 
effect is explored in §3. In §4 we model H11 regions ionized by the composite spectra of starbursts 
of different ages, and investigate the observational consequences of superposed Hn regions (see, 
also, VGS). Since it is the structure of the "outer" zones in an H11 region which most influences 
the icf, in §5 the predictions of radiation-bounded and matter-bounded models are compared. 
Our results are discussed in §6 and our conclusions presented in §7. 
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2. Evolution of Starburst Spectra And The icf 

Cid-Fernandes et al. (1992) explored the time-evolution of the radiation from starbursts of 
different total stellar masses. The models utilize the Maeder stellar evolution models and adopt 
a Salpeter IMF (slope = 1.35) with an upper mass limit of 130 M & . In Figure |l| are shown the 
spectral flux distributions for a starburst of 10 6 M Q at five different times in its evolution: t = 0.0, 
2.5, 3.3, 4.5, 5.4 Myr. To more clearly expose those aspects of these spectra which are key to 
understanding the evolution of the icf, in Figure |2| are shown the time-evolution of various ratios 
of the hydrogen (H°) and helium (He , He + ) ionizing fluxes. 




E (ev) 

Fig. 1. — The spectral distribution (photons s eV ) of the radiation from a 1O 6 M starburst at 
five different epochs in its evolution: t = (solid), 2.5 Myr (dotted), 3.3 Myr (long-dashed), 4.5 
Myr (short-dashed), and 5.4 Myr (dot-dashed). The breaks in the spectra occur at the hydrogen 
and helium ionization edges, indicated by the arrows. 

Initially, due to the presence of massive, hot stars, relatively large amounts of helium ionizing 
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AGE (Myr) 

Fig. 2. — The evolution of the ionizing fluxes as the starburst ages. The solid line traces the 
evolution of hydrogen-ionizing photons {E > 13.6 eV). The dotted line follows the ratio of all He- 
ionizing {E > 24.6 eV) to H-ionizing fluxes, while the dashed line is the same for the He + -ionizing 
photons (E > 54.4 eV) alone. 

photons are present resulting, in general, in icf < 1 (see VGS). As the more massive stars evolve 
and die, the overall flux from the starburst decreases with the relative number of helium ionizing 
photons falling even faster, so that the icf increases and, by t = 2.5 Myr, icf > 1 (VGS). However, 
as the starburst continues to age, stellar winds strip away the atmospheres of the massive stars, 
exposing their hotter interiors. As a result, both the overall flux and the relative numbers of 
helium ionizing photons increase. In these more evolved Hn regions (t ~ 3 Myr) the flux of 
He + -ionizing photons is significant and detectable amounts of doubly-ionized helium should be 
present. As the star cluster continues to evolve, the overall flux of ionizing radiation decreases 
and so, too, does the ratio of helium-ionizing to hydrogen-ionizing photons. From Fig. ^ it can 
be anticipated that while a ^ 5 Myr old Hn region will have an icf > 1, similar to that of a 
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~ 2.5 Myr starburst, detectable amounts of He ++ may be present in the former while absent in 
the latter, perhaps permitting an observational discrimination between the two. 



2.1. Photoionization Models 

Here, as in VGS, the AANGABA photoionization code (Gruenwald & Viegas 1992) is used 
to model the gas whose distribution is taken to be spherically symmetric and homogeneous. 
Since we are specifically interested in the helium - hydrogen icf for low-metallicity Hn regions, 
a metal-poor chemical composition (0.1 solar) is chosen. Initially He/H = 0.083 (Y ~ 0.25) is 
adopted, but later (§3) the consequences for the icf of different helium abundances are explored. 
The radiation intensity is characterized by the number of ionizing photons above the Lyman 
limit, Qh, which is directly related to the mass and the IMF of the stellar cluster. Along with 
the gas density, n, Qh defines the model Hn regions. The popular ionization parameter, U, is 
proportional to Qjj/n, so that a grid of models at fixed density, characterized by different choices 
of U and Rj (the inner radius of the model H n region) , corresponds to the identical grid labelled 
by different values of Qh- 

In VGS the icf (called "ICF" in that paper) was defined by relating the total He/H ratio to 
that of He + /H + alone. As a result, any He ++ was included along with unseen neutral H and He. 
However, lines from the recombination of He ++ are observed in many Hn regions, so that a more 
generally applicable definition, adopted in this paper, is 

He _ . , (He + + He++) 

IF = lcf m • (2) 

For Hn regions with He ++ <C He + , icf ~ ICF, while in the presence of small, but observable 
amounts of He ++ , icf ^ ICF. 

Using AANGABA, grids of model Hn regions ionized by the radiation from starbursts of 
differing masses (i.e., different Qh), were created at each of the cluster ages identified above. In 
Figure ^ the results are shown in the icf - log n plane. As anticipated from Figures [l] and ^, 
icf < 1 for the hard spectra of the zero-age, 3.3, and 4.5 Myr starbursts, while icf > 1 for the 
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Fig. 3. — The icf (see eq. |2|) versus the radiation softness parameter log rj relations for Hn regions 
of all ages (0.0 < t < 5.4 Myr). The solid line is for t = 0, the long-dashed line is for 3.3 Myr, the 
short-dashed line is for 4.5 Myr, and the dot-dashed lines are for 5.4 Myr. The 2.5 Myr model, for 
which log r) ^ 0.8, is off scale. The two curves for the 5.4 Myr starburst are for two choices of the 
helium abundance, He/H = 0.075 and 0.083 (see §||). 

soft spectra of the 2.5 and 5.4 Myr models. Note that the 2.5 Myr starburst is offscale in Fig. ^ 
(where log rj < 0.5) since for the very soft spectrum of this starburst, log r\ ^ 0.8 (see VGS and 
Fig. ||); for the IT data, log r\ < 0.4. For the 5.4 Myr starburst the icf — r\ relation depends on the 
helium abundance as discussed below in §||. For the t = 0, 3.3, 4.5 Myr models, and the 5.4 Myr 
model with He/H = 0.075, as the intensity of the starburst increases, rj decreases and icf — > 1. In 
contrast, for the 5.4 Myr model with He/H = 0.083, as the intensity of the starburst increases, at 
first the icf decreases at nearly constant rj, then r/ increases (to log rj ~ 0.6) at nearly constant 
icf, and finally, as the size of the starburst is increased even further, r\ decreases (to log rj ~ 0.4) 
as icf — > 1. 
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Fig. 4. — The HeII/H/3 ratio versus log rj for Hn regions of differing ages. The starburst ages are 
labelled on the curves. There are two dot-dashed curves for the 5.4 Myr starburst, corresponding 
to the two helium abundances (see Fig. ||). Note that for models of zero age and 2.5 Myr old 
starbursts, the predicted Hell emission is so small as to be off scale. The filled circles are the IT 
data. 

We note that our models compute the line intensities integrated over the whole Hn region 
while the IT data to which we will compare is taken using a 2" x 300" slit. Since the slit is likely 
oriented to include the largest possible fraction of the Hn region, it is reasonable to compare our 
models with their observations. However, since the icf is determined by the outermost zone of 
the Hn region, for those radiation-bounded regions which may be larger than the observing slit, 
the deviation of the icf from unity will be reduced relative to our results. For example, for long 
slit observations of Hn regions smaller than the slit length but larger than the slit width, |1 — icf\ 
can differ by up to 50% from our model predictions. 

In Figure |] the model-predicted HeII(4686)/H/3 flux ratios are shown for starbursts of 
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different ages, as a function of the radiation softness parameter r\. Note that no curves appear 
for the zero-age and 2.5 Myr models because they predict that no detectable He ++ should be 
present (at a level HeII/H/3 ^ 10 -3 ). Also shown are the relevant data from IT who have detected 
He ++ in ~ 30 Hn regions. It is noteworthy that more HeII(4686) is observed than is predicted 
for the young Hii regions, while the opposite appears to be the case for the older Hn regions. 
This strongly suggests that "real" Hn regions may be ionized by a mixture of starbursts of 
different ages or, that what is observed could be the superposition of several H n regions within 
the telescope aperture (see VGS). If so, then the "true" icf must reflect some appropriately 
weighted average of the various icfs shown in Figure [||. We return to this issue in §Q, but note 
here that for the starburst spectrum employed by SJ, it is unlikely that any detectable He ++ will 
be present. Clearly, the presence or absence of He ++ can be a valuable clue to the hardness of 
the ionizing radiation spectrum. Before pursuing these issues, we first consider the alternative 
radiation softness parameter proposed by BFM. 



2.2. An Alternate Radiation Softness Parameter? 

There are advantages and disadvantages to using the Vilchez-Pagel radiation softness 
parameter r\ (see eq. [I]). Since r\ requires the ionic ratios, it may be more convenient to employ a 
related parameter, similar to one introduced by Skillman (1989), which is simply a ratio of line 
ratios, 

/(A3727) /(A6312) 
V V(A4959 + A5007) A /(A6717 + A6731) j ' 1 ' 

From our models we find that for the IT data range of log r] < 0.4, r/ and ry are closely related by, 

log rj = -1.82 + 0.88 log T]. (4) 

In our attempt to compare this with the IT line and ionic ratio data we have uncovered several 
inconsistencies between their line intensity ratios and their ionic ratios, which contribute to a 
spread around a best fit relation from the data of 

log r)'(IT) = -1.85 + 0.76 log r). (5) 
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To facilitate comparison with previous analyses, in this paper we choose to utilize the original 
Vilchez-Pagel parameter rj. 

For the relatively hard spectra corresponding to starbursts of t = 0.0, 3.3, and 4.5 Myr, 
r\ decreases monotonically as the intensity of the starburst, Q#, increases (i.e., r] and Q# are 
anticorrelated) . For these models there is also an anticorrelation between the icf (which is < 1) 
and r\ (see Fig. |||). In contrast, the rj - Q# and r\ - icf relations are not monotonic for the softer 
spectra from starbursts of 2.5 and 5.4 Myr (for which icf > 1). Furthermore, for the latter 
starburst the character of the icf — r\ relation changes with helium abundance. 




log ([OIII]/[OI]) 



Fig. 5. — The radiation softness parameter rj versus the [OIII]/[OI] ratio for Hn regions of all ages 
(the line types are as in Figs. |] and ||). Note that there are two curves (dot-dashed) for the 5.4 Myr 
models, corresponding to He/H = 0.075 (lower of the two) and 0.083 (see §3). The filled circles are 
the IT Hn regions with observed He ++ ; the stars are for those regions without He ++ . 

BFM proposed an alternate radiation softness parameter, the ratio of [OIII]A5007 to 
[OI]A6300. This ratio of observed fluxes is independent of the Hn region metallicity. We have 
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tracked this ratio in our models and in Figure ^ the relationship between the Vilchez-Pagel 
radiation softness parameter, rj, and the BFM parameter [OIII]/[OI], is shown for a series of 
Hn region models ionized by starbursts of different ages. The variable along the curves is the 
intensity of the starburst, proportional to Qh- For all starbursts as Qh increases so, too, does 
the [OIII]/[OI] ratio. Figure || exposes a crucial difference between n and [OIII]/[OI]: whereas all 
spectra, whether "hard" or "soft", span the entire range in [OIII]/[OI], the ranges covered by rj 
are different for hard and soft spectra. While rj may not be the ideal radiation softness parameter, 
the [OIII]/[OI] ratio has virtually no sensitivity to the spectral shape; rather, [OIII]/[OI] provides 
a measure of the intensity of the ionizing radiation field. 

There is another problem if the [OIII]/[OI] ratio is used to provide an estimate of the hardness 
of the radiation spectrum. The [OI]A6300 emission- line is produced in the recombination region, 
so its intensity is very dependent on the optical depth at the Lyman limit, tll- We return to this 
point in §@, but simply remark here that for most models, the [OIII]/[OI] ratio decreases by more 
than 3 orders of magnitude as tll increases from of order unity (matter-bounded) to of order 10 4 
(radiation-bounded). In contrast, although the + and S + fractions are each sensitive to tll, f] 
varies much less with tll since it involves their ratio + /S + . 

The relevant IT data are also shown on Fig. |. Notice that the IT H II regions, with —0.3 ^ 
log 7] < + 0.4, lie between the models for starbursts of all ages, but much below that for a 2.5 
Myr starburst for which log n ^ 0.8. This suggests (see VGS) that starbursts of ~ 2.5 Myr, 
with icf > 1, are not dominant in the H11 regions observed by IT. The icfs for starbursts of 
differing ages are shown as a function of the [OIII]/[OI] ratio in Figure |6[ Once again, notice that 
a measurement of the [OIII]/[OI] ratio does not determine either the size, or the sign, of the icf. 

The predicted HeII/H/3 ratios are shown as a function of [OIII]/[OI] for starburst models of 
different ages in Figure [7|. As was already seen in Fig. ||, no significant Hell is present in the young 
starbursts (t ^5 2.5 Myr), while higher than observed HeII/H/3 ratios are predicted for the older 
Hn regions. 

Although we are in agreement with BFM that the icf — > 1.0 as [OIII]/[OI] increases, we 
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Fig. 6. — The icf versus the [OIII]/[OI] ratio for Hn regions of all ages. The line types are as in 
Figs. |], [|, and [|. The labels on the dot-dashed curves (5.4 Myr starburst) are the values of He/H. 

emphasize that for the [OIII]/[OI] range of the IT data (1.4 < log([OIII]/[OI]) < 2.9), there 
are significant ionization corrections for starbursts of all ages. This is in contrast to the BFM 
suggestion that icf = 1 when [OIII]/[OI] ^ 300. At the same time, the Pagel et al. (1992) proposal 
that icf = 1 for log rj < 0.9 is unsupported by either our results (see Fig. |3|) or those of SJ. For 
the rf and [OIII]/[OI] ranges of the IT data, significant ionization corrections (at the few percent 
level) - most of them with icf < 1 - are to be expected (see, also, SJ). 

3. The Effect Of The Helium Abundance On The Ionization Structure 

Although models of the low metallicity Hn regions targeted for determining the primordial 
helium abundance almost always adopt low heavy element abundances, little attention is usually 
given to the assumed helium abundance. For example, the helium abundance (by number) is 
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Fig. 7. — The HeII/H/3 ratio versus [OIII]/[OI] for Hn regions of differing ages. The starburst ages 
are labelled on the curves. Note that for models of zero age and 2.5 Myr starbursts the predicted 
Hell emission is so small as to be off scale. The filled circles are from IT. 

He/H = 0.10 in Stasinska's (1990) models which were used by Pagel et al. 1992 to analyze their 
observed, low-metallicity Hn regions. This helium abundance, Y « 0.29, is even higher than the 
solar value (Bahcall, Pinsonneault, & Basu 2001). The Armour et al. (1999) and BFM analyses 
use the Orion abundance, He/H = 0.095 (Y » 0.28). VGS set Z = Z /1O and chose He/H = 0.083 
(Y » 0.25). 

Since helium competes with hydrogen for some of the same ionizing photons, the precise 
value of the helium abundance may affect the ionization structure of the H n region (Gruenwald & 
Viegas 2000). As we are tracking small deviations in the derived helium abundance, it is important 
to explore the effect of the adopted helium abundance on the model predictions of the icf, along 
with the other diagnostics such as rj and the [OIII]/[OI] ratio. For starbursts of all ages model 
Hn regions were computed for helium abundances varying from He/H = 0.075 (Y « 0.23) to 0.10 
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(Y ~ 0.29). From these models the following picture emerges. 

The Hn regions formed by massive starbursts (high Qh) are large, with very thin transition 
zones, leading to an icf — ► 1. For the more realistic starbursts (smaller Qh), the He and H° zones 
are larger and the icf deviates from unity, > 1 or < 1 , depending on the hardness of the radiation 
spectrum. With the exception of the 5.4 Myr starburst, there are no significant variations of the 
icf (or of the other diagnostics) with He/H. 
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Fig. 8. — The radial distributions of the H + (solid), He + (dashed), and He ++ (dot-dashed) fractions. 
R is the radius at which the H + fraction has dropped to 10 -4 . The two upper panels are for a 
10 3 Mq starburst and the two lower panels are for a 1O 6 M starburst. In each pair of panels, the 
upper panel is for He/H = 0.075 and the lower one for He/H = 0.100. 

For the softer spectra of the oldest starbursts (^ 5 Myr) there are relatively fewer helium 
ionizing photons (see Fig. ^) and the corresponding Hn regions have some neutral helium inside 
the ionized hydrogen zone so that icf > 1 (see Figs. || & ^). At fixed He/H, as the size of the 
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starburst is reduced (smaller Qh) the transition zone thickens and the relative contribution of 
He increases, leading to an increase in the icf. At the same time, the [OIII]/[OI] ratio (which 
tracks Qh) decreases, leading to an anticorrelation between the icf and [OIII]/[OI]. Now, at 
fixed Qh, as the He/H ratio increases the [OIII]/[OI] ratio, which depends mainly on Qh, is 
essentially unchanged, but the transition zone, where helium has become neutral while hydrogen 
is still ionized, grows in size since there is now more helium competing for the same ionizing 
photons. This effect is shown in Figure || where the radial distributions of the H + , He + , and He ++ 
fractions are shown for two choices of the size of the starburst (M* « 1O 3 M and ~ 1O 6 M0) 
and two choices for the helium abundance (He/H = 0.075, 0.100; Y ~ 0.23, 0.29). For the smaller 
starburst, log([OIII]/[OI]) rs 2, while for the larger one, log([OIII]/[OI]) rs 3 for all values of He/H 
from 0.075 to 0.100. As a result, for a 5.4 Myr starburst, when the helium abundance increases 
the curves shift upwards in the icf - [OIII]/[OI] plane (see Fig. g). 

A much more dramatic change in the character of the curves occurs in the icf - rj plane. For 
He/H < 0.083, as Qh increases both rj and the icf decrease: rj and the icf are positively correlated 
(while rj and [OIII]/[OI] are slightly anticorrelated; see Fig. |5|). In contrast, for He/H ic 0.083 the 
variations of rj and the icf with He/H, and with the size of the starburst, are non-monotonic. For 
these higher helium abundances (Y ^ 0.25), not entirely relevant to our study of low-metallicity 
Hn regions, the structure of the transition zone is very sensitive to the precise helium abundance. 
For example, for the He/H = 0.083 case shown in Fig. [3|, as Qh increases from smaller to larger 
starbursts, the icf first decreases at nearly constant (slightly decreasing) rj. But, as Qh continues 
to increase, and before the icf — ► 1, rj begins to increase at nearly constant icf. Finally, for 
values of Qh larger than that for any observed Hn regions (i.e., [OIII]/[OI] ^> 10 3 ), rj once again 
decreases and so, too, does the icf (icf — > 1). 

4. H II Regions Ionized By Multiple Starbursts 

H II regions ionized by soft spectra will have icf > 1 and relatively high values of rj, while 
those ionized by harder spectra will have icf < 1 and smaller rj values. The IT Hn regions have 
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relatively small rj values (—0.3 < log rj < +0.4; see Fig. ||), confirming the VGS and SJ conclusions 
that by neglecting the icf (i.e., by assuming that the icf = 1) IT have likely overestimated the 
primordial helium abundance. The VGS and SJ results receive further support from Figure |5] 
where the majority of the IT data points are seen to lie between the zero-age model and the 
t = 3.3 — 4.5 Myr starbursts which have harder spectra and even smaller icfs. However, the spread 
in the data highlights the likelihood that at least some Hn regions may be ionized by spectra 
from more than one starburst; there is support for this possibility from the Conti & Vacca (1994) 
and the Conti, Leitherer, & Vacca (1996) observations of substructure within Hn regions. Further 
support comes from the IT detection of doubly- ionized helium in nearly 3/4 of their Hn regions 
at a level above the zero-age prediction, but below the predictions of the 3.3 - 4.5 Myr models (see 
Figs. 4 & 7). Since models of all three of these starbursts predict icf < 1, composite models built 
from these spectra will have intermediate values of rj, and icf < 1. However, what will happen 
if, to these spectra, contributions are added from the 2.5 and 5.4 Myr starbursts? Clearly, the 
resultant Hn regions will have intermediate values of rj, and an icf which could either be < 1 or 
> 1. 

We have explored two types of composite models to see when H n regions ionized by multiple 
starbursts will have sufficiently small values of r? (and of [OIII]/[OI]), but an icf which exceeds 
unity. In VGS we studied composite models of separate Hn regions (ionized by t = and 2.5 
Myr starbursts) which may be included within the telescope aperture, using various line ratios to 
exclude such contamination from the IT data (see VGS for details). Here we extend this analysis 
to combinations of Hn regions ionized by t = and 5.4 Myr starbursts; starbursts of intermediate 
age all have icf < 1 and will only serve to reduce the icf from its t = value. However, before 
addressing such models, first consider H11 regions ionized by the combined spectra of more than 
one starburst. These multiple starburst possibilities are limitless, so we have chosen to investigate 
the most "conservative" options: combining a zero-age starburst with either a 2.5 Myr or a 5.4 
Myr starburst. Including the 3.3 and/or 4.5 Myr starbursts would only further reduce 77 and the 
icf, as well as increase the predicted amount of He ++ (which may already be too high). 
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Despite some similarities between the 2.5 Myr and 5.4 Myr spectra, that of the 2.5 Myr 
starburst is "softer", leading to Hn regions with larger values of rj and of the icf. The key question 
is, what is the resulting icf when the contribution from a zero-age starburst (added to a 2.5 Myr 
starburst) is "large enough" to have reduced r\ to the IT range (log rj < 0.4)? The results, subject 
to the constraint [OIII]/[OI] < 10 3 , are shown in Figure [J It is the [OIII]/[OI] constraint that sets 
the low rj limits and eliminates any composite models with log rj < 0.4 and icf > 1, in agreement 
with our previous conclusions in VGS. 

1.02 | , , , , , 1 , , , , , , , 1 
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Fig. 9. — The icf —log rj relation for a series of Hn regions ionized by the composite spectra for 
zero-age and 2.5 Myr starbursts, subject to the constraint that [OIII]/[OI] < 10 3 . The numbers 
close to the curves are the fractions of the total flux coming from the zero-age starburst. 

The icf versus rj relations for several combinations of the t = 0.0 and t = 5.4 Myr models are 



shown in Figure |10j. Not until the contribution from the aging starburst is some 50 times stronger 
than that from the infant starburst, does the icf exceed unity (and then, only for log rj ^ 0.1). 
However, with such a dominant 5.4 Myr old starburst, significant amounts of doubly-ionized 
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helium may be present. The HeII/H/3 predictions for these same composite models are shown 



in Figure 11. As anticipated, those composite models sufficiently dominated by the 5.4 Myr 
starburst to have icf > 1, will have more Hell than is observed in nearly all of the IT Hn regions. 
Consistent with our earlier results (VGS), as well as those of SJ, it therefore seems unlikely that 
more than a handful of the IT Hii regions can have icf > 1. 




Fig. 10. — The icf— log rj relation for a series of Hn regions ionized by the composite spectra for 
zero-age and 5.4 Myr old starbursts, subject to the constraint that [OIII]/[OI] < 10 3 . For the 5.4 
Myr starburst, He/H = 0.075 is adopted. The numbers close to the curves are the fractions of the 
total flux coming from the zero- age starburst. 

In agreement with SJ, we note from our exploration of these composite models that the 
dividing line between icf < 1 and > 1 corresponds to a ratio of helium-ionizing to hydrogen- 
ionizing fluxes of ~ 0.15. Since our spectra for the individual as well as the composite starbursts 
differ in detail from those of SJ, this result provides supporting evidence that this ratio, and not 
the detailed shape of the ionizing spectra, is key to determining the sign of icf — 1. 
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Fig. 11. — The HeII/H/3— log rj relation for Hn regions with the same He/H, ionized by the same 



composite spectra as in Fig. 10 for the zero-age and 5.4 Myr starbursts. The filled circles are the 
IT data. 



For separate Hn regions, ionized by t = and 5.4 Myr starbursts, which may be included 
within the telescope aperture, the situation is less straightforward. In our exploration of such 
composite models we have found it possible that the combined Hn region can have an icf > 1 
while satifying log rj ^ 0.4 along with the [OIII]/[OI] and HeII/H/3 constraints. Such models 
always have log r\ > 0.25. A glance at Figure [| reveals that only five of the forty-one IT H n regions 
employed in our analysis lie in this range and, of these, only three have observed He ++ . As a 
result, only ~ 3 of the IT Hn regions (< 10%) might have an icf > 1, while still being consistent 
with the observed values of rj, [OIII]/[OI], and HeII/H/3. 
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5. Matter-Bounded Versus Radiation-Bounded 

As mentioned earlier, since the [OI]A6300 emission line is formed in the recombination region, 
its intensity depends on the optical depth at the Lyman limit, tll- As tll increases from 1 
(matter-bounded) to 10 4 (radiation-bounded), [01] increases and the [OIII]/[OI] ratio decreases 
by more than 3 orders of magnitude. Since rj is a ratio of ratios, the effect of the Lyman limit 
optical depth largely cancels and rj is relatively insensitive to tll- 

Images of H n regions show that they are far from being homogeneous, and this suggests 
that different "sectors" could have different values of tll- Indeed, several studies indicate that a 
considerable fraction of the ionizing photons escape from many Hll regions (e.g., Oey & Kennicutt, 
1997; Zurita et al. 2000; Relaho et al. 2001). It is then possible that for the IT-observed Hn regions 
used in determining Yp, parts are matter bounded with a very large - local - [OIII]/[OI] ratio, 
while the remaining radiation bounded sectors have much lower [OIII]/[OI] ratios. Thus, while the 
observed [OIII]/[OI] ratio might be large, suggesting that no ionization correction need be made 
(BFM), this may not necessarily be true. Indeed, as has already been seen, for the IT range in 
[OIII]/[OI], the icf is expected to differ from unity. 

To explore the effect of Hn regions with matter-bounded pieces, several test cases were run 
in which the Hn region was divided into two sectors, a radiation-bounded part with [OIII]/[OI] 
<C 300, and a matter-bounded piece with tll < 10. For model Hn regions with matter-bounded 
sectors having tll ~ 10, ratios of [OIII]/[OI] > 300 can be achieved for covering factors of the 
radiation-bounded part in the range 0.3 to 0.8. As a result, Hn regions might have very large 
[OIII]/[OI] ratios because they are ionized by very massive starbusts (in which case icf — > 1), or 
they may have matter-bounded sections and result from much smaller starbursts. For two-sector 
models, as the contribution from the matter-bounded piece increases, the [OIII]/[OI] ratio 
increases while rj shifts to slightly higher values (compared to the corresponding radiation-bounded 
model). In these cases the icf, although tending towards unity, remains very close to that 
of the radiation-bounded model provided that [OIII]/[OI] < 300. Thus, although partially 
matter-bounded Hn regions might be observed with very high [OIII]/[OI] ratios, suggesting an 
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icf — > 1, the observed values of [OIII]/[OI] and 77 provide some insurance against being lured into 
this trap. Since the [OIII]/[OI] line ratio depends on the number of ionizing photons and on the 
fraction of the solid angle in the outer parts of the H 11 region which are radiation-bounded, it is 
important that observers estimate the total flux in H/3 and the fraction of the flux included in the 
slit. This would permit a better comparison with the models and, by comparing the observed 
Qh for the whole Hn region and the observed [OIII]/[OI] with the model predictions it might be 
possible to distinguish density-bounded H 11 regions from those which are radiation-bounded. 

The discussion here, as well as that in §3, provides a reminder of the importance of a careful 
treatment of the recombination transition zone to accurate, quantitative predictions of the model 
H 11 region icf — 77 relations. In this connection we note that S J assume that the effects of 
radiative transfer can be treated in the "On-The-Spot" approximation (OTS). We have tested this 
assumption and find it to be unjustified. For the spectra of our 0.0, 3.3, and 4.5 Myr models, at 
fixed 77 the OTS underestimates the icf (the icf is too small) by a few percent. In contrast, for 
the 2.5 Myr starburst the OTS overestimates the icf by ^ 0.25. For the 5.4 Myr starburst the 
differences are even more dramatic. For this model, when He/H = 0.083 the icf — 77 relation shown 
in Fig. 3 is shifted upwards (the icf is overestimated), similar to the shift for the 2.5 Myr model, 
but in this case by > 0.3. However, for this same starburst, when He/H = 0.075 the character 
of the icf — 7] relation changes completely from that in Fig. 3 to resemble the one for He/H = 
0.083 and, the icf is shifted upwards to icf > 1.25. Thus, although our results are in excellent 
qualitative agreement with those of SJ, we caution that a careful treatment of radiative transfer is 
crucial for quantitatively accurate, model-predicted Hn region icfs. 



6. Discussion 

In outward appearance, as well as in their internal structure, Hn regions evolve as they age. 
In their infancy H11 regions ionized by young, hot, metal-poor stars have ionized helium present in 
the zone where hydrogen is making the transition from ionized to neutral, resulting in an icf < 1. 
As the starburst ages, the radiation spectrum softens, the icf increases, and in their early youth 
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Hn regions will have icf > 1. So far, no detectable amounts of doubly-ionized helium will be 
present. As Hn regions enter middle age, the radiation spectrum hardens and the icf decreases, 
once again dropping below unity. By now the spectrum has enough sufficiently hard photons that 
detectable amounts of doubly-ionized helium should be present. Finally, as Hn regions enter old 
age, the spectrum softens, the icf increases (icf > 1) one last time, and they fade away. But, 
how does an Hn region show its age? Is there a chronometer? The proposed "radiation softness" 
parameters should be surrogate chronometers but, are they? 

As Figure | reveals, the [OIII]/[OI] ratio is a very poor chronometer. Notice that for 
[OIII]/[OI] ^ 10, the same value of [OIII]/[OI] corresponds to all five spectra studied here. While 
the [OIII]/[OI] ratio provides an excellent measure of the intensity of the ionizing radiation, it 
reveals nothing about the age of the H n region, or the hardness of its spectrum. The situation 
is better for the Vilchez-Pagel parameter rj. Figure |5| suggests that if log rj > 0.7 were observed, 
then the Hn region would likely be dominated by a starburst for which icf > 1. Similarly, it is 
likely that icf < 1 for those Hn regions with log r\ % — 0.1 (see Fig. 5). However, starbursts of all 
ages (except 2.5 Myr) could be responsible for Hn regions with —0.1 % log r\ ^ 0.4. Of these, the 
oldest starburst has an icf > 1, while all the others have an icf < 1. Unfortunately, even in this 
case (low rj), the Vilchez-Pagel radiation softness parameter alone is incapable of distinguishing 
among starbursts of t = 0.0, 3.3, and 4.5 Myr, even though the latter two, older starbursts have 
very different spectra from the youngest one. 

At the same time that Figure || exposes the inadequacy of the [OIII]/[OI] ratio, and the 
deficiencies of rj, it suggests that rj and [OIII]/[OI] in combination might provide a more useful 
gauge of H n region ages by helping to separate the effects of "hard" versus "soft" spectra (which 
determine whether icf < 1 or > 1) from those related to the intensity of the starburst (which 
determines by how much the icf will differ from unity). In addition, for log rj < 0.5, the presence 
or absence of He ++ may help to discriminate between infancy (no observable He ++ at t ~ 0) 
and maturity (detectable He ++ present for 3 ^$ t ^$ 5 Myr); see Fig. ||. For example, Hn regions 
with [OIII]/[OI] ^ 150 and 0.3 < log rj < 0.5 are likely ionized by an old (~ 5.4 Myr) starburst 
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and should have detectable He ++ , while those with the same large values of [OIII]/[OI], but with 
— 0.1 < log T] < 0.2, should be dominated by a very young (t rs 0) starburst and they should have 
no observable He ++ . With reference to Figs. 3, 5, and 6, this approach may be applied to the IT 
data. While many observational constraints are needed for a good photoionization model of a 
specific giant Hn region (SJ), we will account for the range of Hn regions included in the IT data 
set using these three parameters. It is clear that more detailed data for any individual H n region 
may lead to a better model for that specific region, resulting in a slightly different derived value 
for the icf. However, the agreement between our results (below) and those of SJ suggest that 
these differences are likely very small for the bulk of the IT H n regions. 

Of the 41 IT Hn regions in our data set, 36 have log rj < 0.25 (see Fig. 5) and are likely 
ionized by radiation with spectra similar to those of our t = 0.0, 3.3, or 4.5 Myr starbursts. For 
these Hn regions, icf < 1. Of the remaining five Hii regions, only three have He ++ detected; only 
these three are likely to have significant or dominant contributions from a spectrum similar to that 
of our 5.4 Myr starburst for which icf > 1. It seems likely then that for more than 90% (38/41) of 
the IT Hii regions, icf < 1 and Y < Y(IT). This dominance of icf < 1 is not at all surprising when 
considering that an Hii region spends most of its life ionized by a hard radiation spectrum with 
Q(He°)/Q(H°) > 0.15 (see Fig. 2). To estimate the magnitude of the ionization correction it is 
necessary to use the available data for both log r\ and log([OIII]/[OI]) since for the three starbursts 
with hard spectra, the model-predicted icfs vary significantly (see Fig. 3) over the range in n 
covered by the IT data (-0.3 < log r] < 0.25). However, for log([OIII]/[OI]) > 1.7 (see Fig. 5), the 
icfs for these three starburst models are nearly identical and, for log([OIII]/[OI]) ^$ 2.9, they do 
differ from unity (see Fig. 6). As a result of this comparison, we conclude that for ~ 38/41 of the 
IT Hii regions, 0.95 ^ icf ^ 0.99. Since the IT helium abundance determinations range from Y 
« 0.24 to Y ~ 0.26, this suggests a correction to Y P (IT) in the range 0.002 < -AY < 0.010, in 
excellent agreement with the SJ estimate of a ~ 2 — 4 % overestimate in the helium abundance, 
corresponding to 0.005 ^$ —AY ^ 0.010. Approximating our correction as AY ~ —0.006 ± 0.002 
and combining it in quadrature with the IT primordial helium estimate of Yp(IT)= 0.244 ± 0.002, 
leads to a revised estimate of Yp = 0.238 ± 0.003. 



7. Conclusions 



The history of the study of the helium-hydrogen icf extends over several decades (see VGS 
for references to early work) and the different model assumptions may have obscured the fact 
that all such studies reach the same general conclusions. Hn regions ionized by massive, intense 
starbursts will be large and, when radiation bounded, the transition zone from ionized to neutral 
gas is relatively thin compared to the overall size of the H n region. For such H n regions the 
icf — > 1. As the intensity of the starburst diminishes, the relative size of the transition zone 
increases and the icf deviates from unity. Whether the deviation will be icf < 1 or icf > 1 
is determined by the relative hardness of the spectrum of the starburst, not by its intensity. 
Although there will be some quantitative differences among the predicted icfs depending on the 
details of the spectra adopted (e.g., compare VGS, BFM, and SJ), it is clear that the ratio of 
helium-ionizing to hydrogen-ionizing photons is key to determining whether the icf < 1 or > 1. 
Thus, although the spectra adopted here for starbursts of differing ages differ somewhat from the 
corresponding spectra in SJ, we agree with them that there is a "critical" value of Q(He°)/Q(H°) 
~ 0.15 separating Hn regions with icf < 1 from those with icf > 1. An ideal radiation softness 
parameter would provide a unique reflection of this ratio. The Vilchez-Pagel parameter r], although 
imperfect, is a valuable surrogate. However, even when it can be determined that the icf < 1 (or, 
> 1), the amount by which the icf deviates from unity depends on the intensity of the radiation 
spectrum. For more intense (more massive) starbursts, rj decreases, [OIII]/[OI] increases, and the 
icf — > 1. However, the majority of the Hn regions used to infer the primordial helium abundance 
have values of r\ and [OIII]/[OI] which lie in those ranges where the models suggest significant 
(i.e., few percent) deviations of the icf from unity. Note that for an uncorrected helium mass 
fraction in the range Y ss 0.24 - 0.26 (IT), a 3% deviation in the icf from unity corresponds to a 
change in helium abundance of order 0.006, three times as large as the IT error on their inferred 
primordial helium mass fraction. Indeed, in our previous analysis we estimated that the IT value 
of Yp = 0.244 ±0.002 needed to be reduced by approximately 0.003 (VGS). The analysis presented 
here, in agreement with that of SJ, strengthens this conclusion and suggests that the reduction 
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may be even larger. 

The IT primordial helium estimate (Yp = 0.244 ± 0.002) is already somewhat low in 
comparison to the BBN-predicted value based on the baryon density inferred from the deuterium 
observations of O'Meara et al. (2001): for D/H = 3.0 ± 0.4 x 10~ 5 , Q B h 2 = 0.020 ± 0.002, and 
Yp(BBN; D/H) = 0.247 ± 0.001. If, instead, the CMB estimate of the baryon density is used, 
£lsh? = 0.022 ± 0.002 (e.g., Kneller et al. 2001 and references therein), the predicted primordial 
helium abundance is even larger: Yp(BBN; CMB) = 0.248 ± 0.001. Here we have corrected the IT 
estimate for unobserved neutral hydrogen and/or helium, concluding that Yp = 0.238 ±0.003, thus 
increasing the tension between the expected and inferred primordial helium abundances. Almost 
certainly, the discussion here will not be the final word on the primordial abundance of helium. 
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